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We have calculated parity nonconserving 7s — 6d-^^2 amplitude -EpNC in ^23Ra+ using high- 
precision relativistic all-order method where all single and double excitations of the Dirac-Fock 
wave functions are included to all orders of perturbation theory. Detailed study of the uncertainty 
of the parity nonconserving (PNC) amplitude is carried out; additional calculations are performed 
to estimate some of the missing correlation corrections. A systematic study of the parity conserving 
atomic properties, including the calculation of the energies, transition matrix elements, lifetimes, 
hyperfine constants, quadrupole moments of the 6d states, as well as dipole and quadrupole ground 
state polarizabilities, is carried out. The results are compared with other theoretical calculations 
and available experimental values. 

PACS numbers: 31.15.ac, ll.30.Er, 31.15.ap, 31.15.ag 
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I. INTRODUCTION 



There are two separate reasons for parity violation 
studies in an atom: to search for new physics beyond the 
standard model of the electroweak interaction by precise 
evaluation of the weak charge Q.^,, and to probe parity 
violation in the nucleus by evaluating the nuclear anapole 
moment. The atomic-physics tests of the standard model 
that are completed to date were carried out by comparing 
experimental weak charges of atoms Qwi which depend 
on input from atomic theory, with predictions from the 
standard model [l[ . The most precise experimental study 
to date, a 0.35% measurement in Cs was carried out by 
the Boulder group [2] using a Stark interference scheme 
for measuring the ratio of the PNC amplitude E'pnc and 
the vector part of the Stark-induced amplitude /3 for tran- 
sitions between states of the same nominal parity. The 
value of the weak charge in Cs was ultimately found to be 
consistent with the theories of the standard model. How- 
ever, such comparisons provide important constraints on 
its possible extensions. A recent analysis [Sj] of parity- 
violating electron-nucleus scattering measurements com- 
bined with atomic PNC measurements placed tight con- 
straints on the weak neutral-current lepton-quark inter- 
actions at low energy, improving the lower bound on the 
scale of relevant new physics to ^1 TeV. 

Experimental measurements of the spin-dependent 
contribution to the PNC 6s 7s transition in ^'^^Cs led 
to a value of the cesium anapole moment that is accurate 
to about 14% |2i]. The analysis of this experiment, which 
required a calculation of the nuclear spin-dependent PNC 
amplitude, led to constraints on weak nucleon-nucleon 
coupling constants that are inconsistent with constraints 
from deep inelastic scattering and other nuclear experi- 
ments, as pointed out in [3]. Therefore, new experiments 
(and associated theoretical analysis) are needed to re- 
solve the issue. Currently, a microwave experiment to 
measure the spin-dependent PNC amplitude in the 7s 



state of Fr 5] and an isotopic chain experiment in Yb 
are underway. We note that when an experimental 
study is conducted in a single isotope, both theoretical 
and experimental determinations of PNC amplitudes are 
required while the experiments conducted with isotopic 
chains should allow to remove the dependence on the the- 
ory. However, accurate theoretical values for a number of 
atomic properties are useful for this type of experiments 
as well. 

The present work is motivated by the project that was 
recently started at the Accelerator Institute (KVI) of the 
University of Groningen 7] to measure PNC amplitude 
in a single trapped radium ion. Ra+ is a particularly 
good candidate for the PNC study owing to high value 
of the nuclear charge Z and, correspondingly, large ex- 
pected PNC effects. The 7s — 6^3/2 transition in Ra+ 
is of special interest owing to the long life of the 6(^3/2 
state and its sensitivity to both spin independent PNC 
and spin dependent PNC [1]. The 7s — 6(^3/2 transition 
in Ra+ is also being considered for the development of 
optical frequency standards at the same laboratory 
The parity violation experiments are also accompanied 
by a number of measurements of parity-conserving quan- 
tities; as a result we have included a systematical study 
of such properties in this work. 

In summary, we have calculated the PNC amplitude 
for the 7s — 6^3/2 transition in ^^^Ra"*" together with the 
lifetimes of the 7p and 6d states, energy levels for ns, np, 
nd, and nf states, transition matrix elements for a num- 
ber of the El and E2 transitions, quadrupole moments 
of the 6d states, ground state dipole and quadrupole po- 
larizabilities, and magnetic-dipole hyperfine constants A 
for the 7s, 7p, and 6c? states using the relativistic all- 
order method. The all-order method has proved to be 
very reliable for calculating the properties of alkali- metal 
atoms and singly ionized monovalent ions (see, for exam- 
ple, Refs. [10, 11, 12, 13, 14, 15, 16]). The effect of Breit 
interaction on the PNC amplitude is also evaluated. The 
sensitivity of the PNC amplitude to the nuclear radius 
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and varying neutron distribution has been studied. Our 
results are compared with other theoretical values and 
available experimental data. 

II. THEORY 

In this section, we briefly discuss the all-order method 
which has been used to calculate the wave functions and 
the matrix elements necessary to evaluate the observed 
properties. The all-order method relies on including all 
single and double excitations of the core and valence elec- 
trons from the lowest-order wave function: 

ma 7nnab 

+ ^ Pmvalr^av + ^ Pmnvaalaaiaaay]\^y) . (1) 

m^v mna 

Here, is the lowest-order atomic wave function taken 
to be the frozen-core DF wave function of a state v] a|, Qj 
are single-particle creation and annihilation operators, 
Pma and pmv are the single core and valence excitation 
coefRcients, and Pmnab and Pmnva are double core and 
valence excitation coefhcients, respectively. Indices at 
the beginning of the alphabet, a, h, ■ ■ ■ , refer to occupied 
core states, those in the middle of the alphabet m, n, • • • , 
refer to excited states, and index v designates the valence 
orbital. 

To derive equations for the excitation coefhcients, the 
all-order wave function ^ is substituted into the many- 
body Schrodinger equation = E\^v), and terms 
on the left- and right-hand sides are matched, based on 
the number and type of operators they contain. Hamil- 
tonian H = Hq + Vj is taken to be the relativistic no-pair 
Hamiltonian: 

N 

1=1 

Vi = ^"^g^jik ■■ ala]aiak (2) 

ijkl 

where Si are the single-particle energies, : : designate nor- 
mal ordering of the operators with respect to closed core, 
and gijki are the two-body Coulomb matrix elements. 
The all-order equations are solved numerically using a 
finite basis set of single-particle wave functions which 
are linear combinations of B-splines. We have used 70 
basis set B-spline orbitals of order 8 defined on a non- 
linear grid with 500 points within a spherical cavity of 
radius 80 a.u. A large spherical cavity is needed to ac- 
commodate all the valence orbitals required for our cal- 
culation. A sufficiently large number of grid points were 
enclosed within the nucleus to accommodate the influ- 
ence of the nucleus on certain atomic properties such as 
parity- violating matrix elements and hyperfine constants. 

The resulting single-double (SD) excitation coefficients 
are used to calculate matrix elements of various one-body 



operators represented in the second quantization as Z = 

Substituting the expression for the wave function from 
Eq.([T]) in the above equation and simplifying, we get 

_ + + ■ ■ ■ + 

where is the lowest-order DF matrix element and 
Z^°'\ • • • , Z'^') and normalization terms Ni are linear or 
quadratic functions of the single and double excitation 
coefficients [13, [S] • The expression in Eq. ([4]) does not 
depend on the nature of the operator Z, only on its rank 
and parity. Therefore, all matrix elements calculated in 
this work (El, Ml, E2, hyperfine, and PNC matrix ele- 
ments) are calculated using the same general code. 

Corrections to the all-order equations from the 
dominant class of triple excitation terms are also 
evaluated where needed by including the term 
Ipmnrvatalna'lal.avabaal^v) into the SD wave func- 
tion ID) and considering its effect on the energy and 
single valence excitation coefficient equations pertur- 
batively (SDpT approach). Other classes of triple and 
higher excitations are included where needed using 
the scaling procedure by multiplying single excitation 
coefficients p^v by the ratio of the "experimental" and 
corresponding (SD or SDpT) correlation energies p^ . 
The "experimental" correlation energies are determined 
as the difference of the total experimental energy and 
the DF lowest-order values. The calculation of the 
matrix elements is then repeated with the modified 
excitation coefficients. We refer the reader to the review 
[T6| and references therein for the detailed description 
of the all-order method and its extensions. The various 
atomic properties calculated using the all-order method 
described above are discussed in detail in the following 
sections. 



III. PROPERTIES OF RA+ 

A. Energies 

Results of our calculations of energies for a number 
of Ra"*" levels are summarized in Table HI The first 
six columns of Table U give the lowest-order DF ener- 
gies E^^\ the all-order SD energies E^^, the part of 
the third-order energies omitted in the SD calculation 
-E^xtraj first-order Breit contribution B^^\ second-order 
Coulomb-Breit i?^^-* corrections, and Lamb shift contri- 
bution, Ei^s (see Ref. [11] for detail). We take the sum 
of these six contributions to be our final all-order results, 
E^i^ listed in the seventh column of Table [D 
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TABLE I; Contributions to the energies of Ra II; lowest-order (DF) E'-°\ single-double Coulomb all-order correlation energy 
E^^ , third-order terms not included in the SD value -Eg^Lai first-order Breit and second-order Coulomb-Breit corrections B^"\ 
and Lamb shift E ls- The total energies -E^^^t ^^re compared with experimental energies i?oxpt p^l. [Tsl]. 5E = -B^?t - ^'cxpt- Our 
predicted energy values are listed for the 9pi/2 and lOpj energy levels in separate rows. Units: cm"'^. 



nlj 






extra 




5(2) 


-Els 






5E^° 


7si/2 


-75898 


-6692 


1152 


147 


-250 


33 


-81508 


-81842 


334 


6d3/2 


-62356 


-8042 


1152 


155 


-398 





-69488 


-69758 


270 


64/2 


-61592 


-7034 


926 


114 


-360 





-67947 


-68099 


152 


7pi/2 


-56878 


-4027 


587 


102 


-109 





-60326 


-60491 


165 


7fl3/2 


-52906 


-3020 


433 


63 


-90 





-55519 


-55633 


114 


8Sl/2 


-36860 


-1745 


316 


46 


-74 


7 


-38311 


-38437 


126 


7d3/2 


-31575 


-1590 


245 


39 


-92 





-32973 


-33098 


125 


74/2 


-31204 


-1456 


204 


29 


-84 





-32509 


-32602 


93 


5/5/2 


-28660 


-4438 


371 


11 


-63 





-32780 


-32854 


74 


5/7/2 


-28705 


-4159 


353 


8 


-61 





-32564 


-32570 


6 


8pi/2 


-30053 


-1298 


201 


39 


-42 





-31152 


-31236 


84 


8p3/2 


-28502 


-1034 


156 


25 


-36 





-29391 


-29450 


59 


9Sl/2 


-22004 


-741 


136 


21 


-33 


2 


-22618 


-22677 


59 


9pi/2 


-18748 


-605 


96 


20 


-21 





-19259 


















-19305" 






9p3/2 


-17975 


-495 


76 


13 


-18 





-18399 


-18432 


33 


84/2 


-19451 


-683 


105 


18 


-40 





-20051 


-20107 


56 


84/2 


-19261 


-634 


90 


13 


-37 





-19829 


-19868 


39 


lOSl/2 


-14651 


-388 


72 


11 


-18 


1 


-14972 


-15004 


32 


lOpi/2 


-12838 


-335 


53 


11 


-11 





-13120 


















-13144" 






10p3/2 


-12397 


-278 


43 


7 


-10 





-12635 


















-12653" 






94/2 


-13226 


-366 


56 


10 


-22 





-13548 


-13578 


30 


94/2 


-13118 


-342 


49 


7 


-20 





-13424 


-13447 


23 


104/2 


-9587 


-221 


34 


6 


-13 





-9780 






104/2 


-9519 


-207 


30 


4 


-12 





-9704 







" Our predicted values. 



TABLE II: Comparison of the excitation energies important 
to the calculation of the 7s — 64/2 PNC amplitude. All results 
are in cm~^. 



Transition 


Present 


Ref. [20] 


Ref. [21j 


Expt. 


7s - 7pi/2 


21182 


21279 


21509 


21351 


7s - 7p3/2 


25989 


26226 


26440 


26209 


64/2 - 7pi/2 


9162 


9468 


9734 


9267 


64/2 - 7p3/2 


13969 


14415 


14665 


14125 



The column labeled 5E^^ in Table U elves differences 
between our ab initio results and the experimental val- 
ues [13, [3 ■ The SD results are in good agreement with 
the experimental values taking into account very large 
size of the high-order correlation corrections. We predict 
the energies of the 9j5i/2, lOpi/2, and IOP3/2 levels using 
our theoretical results and differences between our and 
experimental values for the known np levels. The pre- 
dicted values are listed in Table |T] and are expected to be 
accurate to a few cm^^. 

We compare our results for the excitation energies im- 
portant to the calculation of the 7s — 6^3/2 PNC am- 
plitude with other theoretical calculations and experi- 



ment [H] in Table HJ The calculations in both Ref. ^ 
and Ref. [2l[ use high-precision all-order methods, but 
represent very different approaches. The calculations 
in Ref. f20| are performed using the correlation poten- 
tial method. The results of Ref. [2l| are obtained us- 
ing coupled-cluster method including single, double, and 
partial triple excitations. The results of Ref. [2(3| are in 
better agreement with experiment for the 7s — 7p transi- 
tions and the results from the present work are in better 
agreement with experiment for the 6(^3/2 — 7p transitions. 
Large discrepancies of the coupled-cluster results from 
Ref. [m for the 6d — Ip transitions with experiment are 
somewhat surprising and may indicate insufficient num- 
ber of higher partial wave functions in the basis set. In 
our calculations, all partial wave up to Imax = 6 are 
explicitly included in all calculations and extrapolation 
for higher number of partial waves is carried out for the 
dominant second-order correlation energy contribution. 



B. Electric-dipole matrix elements 

We calculate all allowed reduced electric-dipole matrix 
elements between ns, np, and nid states where n = 1— 10 
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TABLE III: Comparison of the present results for the absolute 
values of the electric-dipole reduced matrix elements in Ra II 
with other theoretical calculations. All results are in atomic 
units. The lowest-order DF values are listed in the column la- 
beled "DF" to illustrate the size of the correlation correction. 
Negative sign of the DF value for the 8pi/2 — 7si/2 transition 
indicates that the lowest-order value is of the opposite sign 
with the final result. 



Transition 


DF 


Present 


Ref.f20] 


Rcf.[9] 


Ref.[21] 


7jDi /o — 7Si /o 

' L j A ''-'1/iJ 


3.877 


3.254 


3.224 


3.28 


3.31 


7pi/2 - 8Si/2 


2.637 


2.517 


2.534 






7pi/2 - 9Si/2 


0.716 


0.702 


0.708 






7pi/2 - 6^3/2 


4.446 


3.566 


3.550 


3.64 


3.68 


7Pl/2 - 7^3/2 


4.527 


4.290 


4.358 






7pi/2 - 8d3/2 


1.584 


1.445 


1.432 






7jDq /9 — 7Sl 10 


5.339 


4.511 


4.477 


4.54 


4.58 


7P3/2 - 8Si/2 


4.810 


4.644 


4.663 






7p3/2 - 9Si/2 


1.078 


1.035 


1.036 






7p3/2 - 6d3/2 


1.881 


1.512 


1.504 


1.54 


1.56 


7p3/2 - 7d3/2 


2.488 


2.384 


2.407 






7P3/2 - 8d3/2 


0.733 


0.652 


0.641 






7P3/2 - 6^5/2 


5.862 


4.823 


4.816 


4.92 




7P3/2 - 7d5/2 


7.249 


6.921 


6.995 






7P3/2 - 8d5/2 


2.227 


2.011 


1.954 






8pi 10 — 7si 10 


-0.125 


0.047 


0.088 


0.04 




Spi/2 - 8Si/2 


7.371 


6.949 


6.959 






8j)i /o — 9Si /o 


5.227 


5.012 


5.035 






8pi/2 - 6^3/2 


0.105 


0.049 


0.013 


0.07 




8pi/2 - 7d3/2 


10.21 


9.553 


9.540 






8pi/2 - 8d3/2 


7.184 


7.010 


7.104 






8P3/2 - 7Si/2 


0.625 


0.395 


0.339 


0.50 




8P3/2 - 8Si/2 


9.880 


9.294 


9.320 






8p3/2 - 9Si/2 


9.244 


9.022 


9.036 






8P3/2 - 6(i3/2 


0.168 


0.144 


0.127 


0.15 




8P3/2 - 7^3/2 


4.331 


4.035 


4.028 






8P3/2 — 8d3/2 


4.047 


4.002 


4.034 






8P3/2 - 6^5/2 


0.462 


0.378 


0.347 


0.40 




8P3/2 - 7d5/2 


13.37 


12.55 


12.53 






8P3/2 - 8d5/2 


11.68 


11.49 


11.58 







suits by the appropriate angular factors for the pur pose 
of comparison. The calculations of the Refs. |21|| are 
carried out using the coupled-cluster method. 

We have also listed the lowest-order DF values in the 
first column of the table to illustrate the size of the cor- 
relation corrections for various transitions. Negative sign 
of the DF value for the 8pi/2 — 7si/2 transition indicates 
that the lowest-order values is of the opposite sign with 
the final result. The correlation corrections for the pri- 
mary 7s — 7p and 7p — Gd transitions are quite large, 
18-25%. The correlation corrections for the remaining 
strong transitions are generally smaller, 2-10%. All theo- 
retical values are in good agreement for these transitions. 
Our values for 7s — 7p and 7p — 6d are in better agree- 
ment with results of Ref. f29| than those of Refs. 'o', zTl. 
The agreement is generally poorer for the transitions with 
small values of the matrix elements as expected owing to 
very large size of the correlation corrections. Since differ- 
ent methods omit or include somewhat different classes 
of the high-order corrections, discrepancies are expected 
when such corrections are large. The issue of the very 
small matrix elements, such as Bp — 7s, is also discussed 
in Ref. [20]. 

C. Polarizabilities 

We calculate the static dipole and quadrupole polar- 
izabilities of the Ra+ ion in its ground 7s state. The 
static polarizability is calculated as the sum of three 
terms representing contributions from the ionic core Ofc, 
a small counteracting term to compensate for the excita- 
tions from the core states to the valence state a^c, and 
valence polarizability a^,: 

a = etc + a«c + a^- (5) 
1. Dipole polarizability 



and ni = 6 — 10 using the method described above. The 
subset of these matrix elements is compared with the cor- 
relation potential calculations of Ref. [20*1 and coupled- 
cluster calculations of Refs. [1, [2l| in Table Hill Absolute 
values of the reduced matrix elements in atomic units are 
listed in the table. All present values with the exception 
of the 7pi/2 - 8s, 7p3/2 - 8s, 8pi/2 - 7s and 8^3/2 - 7s 
transitions are ah initio SD values. For these four transi- 
tions, we used scaling procedure described above to pro- 
vide recommended values as we expect the scaled values 
to be more accurate based on Cs "best set" data Ref. '23'] . 
The calculations of Ref. [13] are carried out using fitted 
Brueckner orbitals (i.e. include semi-empirical correction 
to the correlation operator) and include core polariza- 
tion, structure radiation, and normalization corrections. 
We note that Ref. [2^ quotes radial integrals rather than 
reduced matrix elements, so we have multiplied their re- 



The valence polarizability contributes over 90% of the 
total value of the electric-dipole polarizability and is cal- 
culated using sum-over-states approach: 

^^(EX) = -y (" K^^ll-Pll^py^)!' + Wsm\np^i2)? \ 

(6) 

The sum over n in Eq. ^ converges extremely fast. In 
fact, the first term with n — 7 contributes 99.8% of the 
total value. As a result, we calculate the first few terms 
(with n = 7 ~ 10) using our all-order matrix elements 
from Table IIIII and experimental energies [TtI . [18| where 
available. The remainder a^*^'' is calculated in the DF ap- 
proximation without loss of accuracy. The ionic core con- 
tribution Q!c and term aye are calculated in the random- 
phase approximation (RPA). The RPA core value is ex- 
pected to be accurate to better than 5% (see Ref. [23] 
and references therein). All contributions to the dipole 
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TABLE IV: Contributions to the ground state dipole polar- 
izability of Ra+. The contributions from the (7 — 10)p states 
are given separately. Our result is compared with calculation 
from Ref. AH results are in a.u. 



Contribution 


CtEl 


7pi/2 - 7s 


36.29 


7P3/2 - 7s 


56.79 


8pi/2 - 7s 


0.00 


8P3/2 - 7s 


0.23 


(9 - 10)p - 7s 


0.04 


„ main 
Q„ 


93.35 


etc 


13.74 


„ tail 

a„ 


0.11 


CX-vc 


-0.98 


Total 


106.22 


Theory [9] 


106.12 



polarizability are listed in Table IIVI The contributions 
from n = 7 — 10 are given together as a™^^^. 

The value of the ground state Ba"*" polarizability cal- 
culated by the same approach l 24l| is in near perfect 
agreement with the experiment [2^ (to 0.2%). Moreover, 
the theoretical SD 6p lifetimes in Ba"*" are also in excel- 
lent agreement with experimental values [23|. We note 
that lifetime experiments are conducted enti rely differ- 
ently from the polarizability measurement of [2^ . There 
are two differences between the Ba"*" and Ra"*" dipole po- 
larizability calculations: increased ionic core contribution 
and increased size of the correlation corrections. The 
core contribution increases from 8% in Ba''" to 13% in 
Ra+ and the correlation correction contribution to the 
7s — 7p matrix elements increases by about 3% (from 
16.6% to 19.1% for the 7s - 7pi/2 transition). Neither 
of these changes is expected to significantly decrease the 
accuracy of the Ra^ ground state dipole polarizability 
in comparison with the Ba"*" one. Therefore, we expect 
our value to be accurate to better than 1%. Our result 
is in agreement with the coupled-cluster calculation of 
Ref. 0. 

2. Quadrupole polarizability 

The valence part of the quadrupole polarizability is 
calculated using the sum-over-states approach as: 

a,„(£2) -^y( l(7^IIQII"^3/2)P ^ m\Q\\nd,/,)\A 

. .. (7) 

All contributions to the quadrupole polarizability are 
listed in Table |Vl The correlation correction to the E2 
matrix elements is dominated by a single term among 
twenty terms in the numerator of Eq. ([4]). As described 
in detail in Ref. [IJ, additional omitted correlation cor- 
rection to this term may be estimated by the scaling pro- 
cedure described above. The scaling modifies the SD re- 



TABLE V: Contributions to the ground state quadrupole po- 
larizability and the E2 reduced matrix elements of Ra"'" in a.u. 
The comparison of our result with other theoretical calcula- 
tion 9] is also presented. 



Contribution 


E2 


QB2 


6^3/2 - 7s 


14.74(15) 


789(13) 


6^5/2 - 7s 


18.86(17) 


1136(16) 


7d3/2 - 7s 


14.21(30) 


182(3) 


7^5/2 - 7s 


16.49(38) 


243(4) 


8^3/2 - 7s 


5.63(4) 


22.6(2) 


8^5/2 - 7s 


6.79(6) 


32.6(2) 


9^3/2 - 7s 


3.30(3) 


7.0(1) 


9^5/2 — 7s 


4.03(3) 


10.4(1) 


10d3/2 - 7s 


2.27(3) 


3.1 


104/2 - 7s 


2.79(3) 


4.7 


„ main 




2430(21) 


„ tail 
Ov 




35(10) 


Oc 




68(12) 


Total 




2533(26) 


Theory [9] 




2547.5 



suits by 0.7 to 2.3% depending on the transition. We 
have also carried out the ab initio all-order calculation 
with inclusion of the triple valence excitation coefficients 
as described in the Section [IT] (SDpT approach). The 
scaling procedure was repeated starting from the SDpT 
approximation for the dominant 7s — 6^3/2 and 7s — 6(^5/2 
transitions. These additional calculations allow us to di- 
rectly evaluate the uncertainty in our calculations since 
they produce different evaluations of the omitted corre- 
lation correction. We take the uncertainty in the calcu- 
lation of the 7s — 6(^3/2 and 7s — 6^5/2 to be the maxi- 
mum of the difference of out final SD scaled results with 
ab initio and scaled SDpT data. We note that SD ap- 
proach generally underestimates the correlation energy 
and SDpT approach generally overestimates the correla- 
tion energy used in the scaling procedure. The scaled 
SD and SDpT results are rather close, further confirm- 
ing the validity of this procedure and of our uncertainty 
estimate. Therefore, we take the uncertainty of the re- 
maining transitions to be the difference of the final SD 
scaled and ab initio SDpT values. The resulting final ma- 
trix elements and their uncertainties are listed in Table IVl 
in column labeled "E2" . The relative uncertainty of the 
corresponding polarizability values is twice the relative 
uncertainty of the matrix elements since we assume the 
experimental energies be accurate to all figures quoted. 
The sum over n converges far slower than in the case of 
the dipole polarizability so calculating a first few terms 
to high precision is essential to obtain an accurate final 
value. The tail contribution, while small, is significant 
and has to be treated with care. We estimated that DF 
value for the main {n — 6 — 10) term is larger than our 
final all-order result by 22%. Therefore, we decrease the 
DF tail of 45 a.u. by 22% and take the difference of the 
DF tail and the final adjusted value to be its uncertainty. 
The core contribution is calculated in the RPA approxi- 
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TABLE VI: Contributions to the lifetimes of the lpi/2 and 
7p3/2 states. The transitions rates A are given in 10^ s^^ and 
the lifetimes are given in ns. 



^(7pi/2 - 

-r(7pi/2) 



7pi/2 



-7s) 
- 6(^3/2) 



104.4 
10.3 
114.7 
8.72 ns 



7p3/2 



^(7p3/2 - 7s) 

^(7p3/2 - 60(3/2) 
^(7p3/2 - 6d5/2) 

■^(7p3/2) 



185.5 
3.3 
22.8 
211.6 
4.73 ns 



mation; we take the difference between DF and RPA core 
values to be the uncertainty of the core contribution. Our 
final value is in agreement with the result of Ref. @ . 



D. Lifetimes of the 7p and 6d states 

The lifetimes r of the 7p and Qd states in Ra+ are cal- 
culated as the inverse of the sum of the transition proba- 
bilities A. The 7p states decay via strong electric-dipole 
transitions. Total of five El transitions contribute to the 
lifetimes of these two states: lpi/2 — 7s, 7pi/2 — 6^3/2, 
7p3/2 - 7s, 7p3/2 - 6(i3/2, and 7^3/2 - 6^5/2. The electric- 
dipole transition rates are calculated using formula 



\E1 



2.02613 X 10^^ \{i\\D\\f)\ 



A3 



2j, 



1 



(8) 



where A is the wavelength of the transition in A and 
is the electric-dipole reduced matrix element 
in atomic units. We use the experimental wavelength 
[itI ITsI I and our all-order matrix elements listed in Ta- 
ble |TTT] when evaluating the transition rates. The results 
are summarized in Table IVTl We find that while the con- 
tributions of the 7s — 7p transitions to the 7p lifetimes 
are dominant, the contributions of the 7p — 6d transitions 
are significant (over 10%). 

Only one transition, 6(^3/2 — 7s, has to be considered for 
the calculation of the 6c?3/2 lifetime. The corresponding 
transition rate is calculated as 



aE2 _ 



1.11995 X 10^ 



A^ 



2J^ 



1 



(9) 



where A is the wavelength of the transition in A and 
(illQII/) is the electric-quadrupole reduced matrix ele- 
ment in atomic units. 

Two transitions have to be considered in the calcula- 
tion of the 60(5/2 lifetime: E2 60(5/2 — 7s transition and 
Ml 6(^5/2 — 6(^3/2 transition. The Ml transition rate is 
calculated as 



AMI 



2.69735 X 1013 \{i\\Ml\\f)\'^ 
A3 2j, + 1 



(10) 



We use the experimental wavelengths [17|, ll8| and our 
all-order matrix elements listed in Table |V] when evalu- 
ating the E2 transition rates. Our result for the reduced 



TABLE VIL Lifetimes of the 6^3/2 and 66(5/2 states of Ra"*" 
in seconds. Comparison of our results with other theoretical 
calculations is presented. 



Term 



r(6c(3/2) 



Present 
Theory [9] 
Theory [20] 



0.638(10) 

0.627(4) 

0.641 



r(6c(5 



/2J 



0.303(4) 
0.297(4) 
0.302 



Ml ~ 60(3/2 matrix element is 1.55 a.u. The E2 and 
Ml transition rates contributing to the 60(5/2 lifetime are 
3.255 s~i and 0.049 s"^. We verified that the contribu- 
tion of the 60(5/2 — 6^3/2 E2 transition is negligible. 

Our results for the 60^3/2 and 6^5/2 lifetimes are pre- 
sented in Table IVIII together with other theoretical val- 
ues. Our values for the lifetimes of the 6o( states are in 
better agreement with those published by Dzuba et al. 
[20| than with the results of Sahoo et al. 9]; however, 
the discrepancies with Ref. 9] are small. We also list 
the uncertainties of our values in Table IVTIl The relative 
uncertainties in our values of the 6d lifetimes are twice 
the relative uncertainties in the values of the E2 matrix 
elements listed in Table IVl We note that the estimated 
uncertainties quoted in Ref. are obtained by carrying 
out calculations with different bases; i.e. they are numer- 
ical uncertainties resulting from the particular choice of 
the basis set and do not include estimation of the miss- 
ing correlation effects. In our calculations, the basis set 
is complete (70 splines for each partial wave) and increas- 
ing its size does not change the result. Our uncertainties 
include estimation of the terms beyond triple contribu- 
tions as described above as well as uncertainty owing to 
truncation of the partial waves above I > 6. Therefore, 
while our uncertainty is higher for 6^3/2 state than the 
one quoted in Ref. [9], it represents an attempt to pro- 
vide an actual boundary for the recommended value of 
this lifetime. 



E. Quadrupole moments of the 6d states 

We also calculated the values of the quadrupole mo- 
ments of the 60(3/2 and 6^5/2 states since these proper- 
ties are of interest to the investigation of possible use of 
Ra"*" for the development of optical frequency standard 
0. The quadrupole moment 0(7 J) can be expressed via 
the reduced matrix element of the quadrupole operator 
Q as 



e(7J) = 



(2J)! 



v/(2J-2)!(2J + 3)! 



(vI/(7J)||Q|lM.(7J)) 



(11) 

The calculation follows that of the E2 matrix elements. 
As in the case of the E2 7s — nd matrix elements, a single 
correlation correction term is dominant, and the omit- 
ted correlation contributions may be estimated via the 
scaling procedure. We have conducted four different cal- 
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TABLE VIII: Quadrupole moments of the 6d3/2 and 6d5/2 
states in Ra""" in a.u. 



State 


SD 


SDpT 




SDpT,, 


Final 


Ref. [9] 


6^3/2 
6^5/2 


2.814 
4.311 


2.868 
4.380 


2.839 
4.342 


2.829 
4.329 


2.84(3) 
4.34(4) 


2.90(2) 
4.45(9) 



TABLE IX: Magnetic-dipole hyperfine constants A(MHz) for 
the 7s, 7pi/2, 7p3/2, 6^3/2, and &d^/2 states in '^■^^Ra"'" calcu- 
lated using SD and SDpT all-order approaches. Lowest-order 
(DF) values are also listed to illustrate the size of the correla- 
tion corrections. The present values are compared with other 
theoretical 0, [2l|] and experimental values from Refs. p7l.[2^. 



State 


DF 


SD 


SDpT 


Ref. [21] 


Ref. [9] 


Expt. 


7s 


2614 


3577 


3450 


3557 


3567 


3404(2) 


6^3/2 


52.92 


81.51 


79.56 


79.80 


77.08 




6^5/2 


19.24 


-23.98 


-24.08 




-23.90 




7Pl/2 


444.5 


699.5 


671.5 


671.0 


666.9 


667(2) 


7p3/2 


33.91 


56.62 


54.40 


56.53 


56.75 


56.5(8) 



rounded off value /i/ = 0.271(2) /ia? from [SPl was used 
in Ref. [9]. The values for A/gj were quoted in Ref. [2l| . 
so we multiplied their values by 0.18067 for comparison. 
The differences between our results and experimental val- 
ues are 1.3%, 0.7%, and 4% for 7s, 7pi/2, and 7p3/2 
states, respectively. We note that the uncertainty in the 
value of the nuclear magnetic moment is 0.7%. Larger 
difference of the ^(7p3/2) SDpT value with the experi- 
ment is similar to that one in Cs , where the difference 
of the SDpT value for the 6^3/2 magnetic-dipole hyper- 
fine constant with experiment is 3.5%. Interestingly, the 
Cs SDpT values are below the experimental ones while 
the Ra+ SDpT results are above the experimental val- 
ues. This can be explained by the uncertainty in the 
treatment of the finite size correction, uncertainty in the 
value of Ra nuclear magnetic moment, and the difference 
in the size and distribution of the correlation corrections 
in Cs and Ra+. 



IV. PARITY NONCONSERVATION 



culations: ah initio SD and SDpT, and scaled SD and 
SDpT ones to evaluate the uncertainty in the final val- 
ues. The results are summarized in Table IVIIII The 
correlation correction to the quadrupole moments is on 
the order of 20%. Our values are compared with coupled- 
cluster calculation of Ref. [9| . Our results are lower than 
that of Ref. This issue has been discussed in detail 
in Ref. where we have demonstrated that CCSD(T) 
method may overestimate quadrupole moments by a few 
percent owing to the cancellation of various terms. Omis- 
sion of orbitals with I > 4 from the basis set may also 
lead to higher values. 



F. Magnetic-dipole hyperfine constants 

Our results for the magnetic-dipole hyperfine constants 
A(MHz) in ^^^Ra+ are compared with theory [i,[2H and 
experiment in Table IIXI The gyromagnetic ratio 

gi for ^^^Ra is taken to be gi = 0.1803 and corresponds 
to the value fij = 0.2705(19) fiN from Ref. We note 
that the magnetic moment of ^^■^Ra have not been di- 
rectly measured but recalculated from measurements of 
^^■^Ra and ^^^Ra nuclear magnetic moments in Ref. (29|. 
The magnetization distribution is modeled by a Fermi 
distribution with the same parameters as our charge dis- 
tribution {c~ 6.862 fm and 10%-90% thickness param- 
eter is taken to be t—2.3 fm). The lowest-order values 
are also listed to demonstrate the size of the correlation 
corrections for various states. The triple contributions 
are important for the hyperfine constants and are par- 
tially included as described in Section |lll These values 
are listed in column labeled "SDpT" . The SD values are 
also listed for comparison in column labeled "SD" . 

The value gi = 0.18067 that corresponds to the 



Nuclear-spin independent PNC effects in atoms are 
caused by the exchange of a virtual Zq boson between 
an electron of the atom and a quark in the nucleus, 
or between two atomic electrons f3d\. The second ef- 
fect is extremely small and will not be considered in this 
work. The dominant PNC interaction between an atomic 
electron and the nucleus is described by a Hamiltonian 
AsVn, which is the product of axial- vector electron cur- 
rent Ae and vector nucleon current Vn- The PNC inter- 
action leads to a non-zero amplitude for transitions oth- 
erwise forbidden by the parity selection rule, such as the 
6^3/2" 7s transition in singly ionized radium. Combining 
experimental measurements and theoretical calculations 
of the PNC amplitude permits one to infer the value of 
the weak charge Qw for precise atomic-physics tests of 
the standard model. 

The 7s— 6^3/2 PNC amplitude in Ra+ can be evaluated 
as a sum over states: 



PNC 



E 



(6^3/2 |i:'|npi/2)(npi/2|gpNc|7s) 

Ets — Enp^,^ 



E°° (6^3/2 1 gpNC |»P3/2) ("-P3/2 1 ^1 7s) 
^6(^3/2 



where D is the dipole transition operator. The values of 
rrLj are customary taken to be rrij — 1/2 for all states. 
The PNC Hamiltonian -ffpNC is given by 



^PNC ^^QwlbP{r) 



(13) 



where is the universal Fermi coupling constant, Qw 
is the weak charge and 75 is the Dirac matrix associated 
with pseudoscalars. The quantity p{r) is a nuclear den- 
sity function, which is approximately the neutron density. 
In our calculations, we model p(r) by the charge form 
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TABLE X: Contributions to the Ep^'g in ^^^Ra^ in units of 10-^H\e\ao{-Qw /N). D and i/pNC are dipole and PNC matrix 
elements, respectively. Reduced electric-dipole matrix elements are listed for consistency with previous tables; they need to 
be multiplied by l/\/6 to obtain relevant values of {i|_D|j) (rrij — 1/2 for all states). All values are in a.u. Our results are 
compared with calculations of Ref . 20] . 



71 


/ / O 1 D\\T},T)^ /o\ 




Er, — Er,r, 




Ref. [20] 


7 


3.566 


-2.665 


-0.0973 


39.882 


40.69 


8 


0.049 


-1.590 


-0.2306 


0.137 


0.11 


9 


0.017 


-1.124 


-0.2849 


0.027 


0.02 


10 


0.008 


-0.841 


-0.3130 


0.009 




n 


(6d3/2|-ffpNcjnp3/2) 


{np3/2\\D\\7s) 




-BpNC 


Ref. [20] 


7 


-0.047 


-4.551 


-0.0644 


-1.348 


-2.33 


8 


-0.040 


-0.405 


-0.1837 


-0.036 


-0.05 


9 


-0.032 


-0.140 


-0.2339 


-0.008 


-0.01 


10 


-0.026 


-0.069 


-0.2602 


-0.003 





factor, which is taken to be a Fermi distribution with 
50% radius Cpnc = Ccharge =6.8617 fm [H and 10%-90% 
thickness parameter t = 2.3 fm for ^^'^Ra"'" , i.e. we take 
p{r) to be the same distribution as the charge distribu- 
tion used our entire all-order calculation of the Ra"*" wave 
functions and corresponding properties. We also investi- 
gate how the PNC amplitude vary with changes in both 

CpNC and Cchargc- 

The sum over n in Eq. (|12|) converges very fast in our 
case, and only first few terms need to be calculated ac- 
curately. Therefore, we divide our calculation of -EpNC 
into three parts: a main term E^^™ that consists of the 
sum over states with n = 7 — 10, a tail -EpNC which is the 
sum over states with n — 11, . . . , cxd, and the contribution 
-^PNC from autoionizing states given by the terms with 
n = 2 — 6. The calculation of the main term is illustrated 
in Table |Xl where we list the "best set" of the dipole 
and PNC matrix elements used in our calculation as well 
as relevant energy differences. The final electric-dipole 
matrix elements are taken to be ah initio single-double 
all-order results (following the comparison of the simi- 
lar Cs and Ba+ results with experiment [111, H^). Re- 
duced electric-dipole matrix elements are listed for con- 
sistency with previous tables; they need to be multiplied 
by l/y^6) to obtain relevant values of {i\D\j) {rrij = 1/2 
for all states). The final PNC matrix elements for the 
6^3/2 ~ 7pi/2 and 6(^3/2 — 8pi/2 transitions are taken to 
be SD all-order scaled values since the contribution that 
can be accounted for by scaling is the dominant one for 
these cases; remaining PNC matrix elements are taken 
to be ab initio SD values. Experimental energies are 
used where they are available, our predicted energy val- 
ues from TableUare used for the 9pi/2, lOpi/2, and lOps/Q 
levels. Our results are compared with results of Ref. |2(j | 
calculated using the correlation potential method. The 
main part of the PNC amplitude is overwhelmingly dom- 
inated by a single term listed in the first row of Table 1x1 
Our result for this term slightly differs from the calcula- 
tion in Ref. [13 (by 2.2% ). However, the Ref. does 
not list the Ra+ isotope for which the calculation has 
been conducted. Since the value of the PNC amplitude 



is multiplied by the neutron number in the present com- 
monly accepted units of 10~^^i|e|ao(— Qi^/iV), the differ- 
ence between values for the PNC amplitudes for ^^^Ra+ 
and ^^^Ra"*" is 2% just owing to 138/135 neutron number 
ratio. Therefore, the difference may be either explained 
by the simple isotope rescaling, difference in the choice 
of the nuclear density function parameters, or differences 
in the treatment of the correlation correction. The only 
significant discrepancy between our calculation and that 
of Ref. 20] is in the other term with n = 7 (-1.35 vs. 
-2.33). This difference has to result from the differences 
in the treatment of the correlation correction since this 
entire value comes from the correlation effects. Taking 
into account that the DF value for this term is consis- 
tent with zero and random-phase approximation (RPA) 
result, -4.08, is larger than the all-order value by nearly 
a factor of 3, such discrepancy is not very surprising. 

To provide some estimate of the uncertainty in the 
calculation of the main term, we conduct the "scatter" 
analysis of the data following the calculation of the Cs 
PNC amplitude [s^]- In such analysis, sets of data for 
dipole matrix elements, PNC matrix elements, and ener- 
gies are varied (i.e., taken to be SD, SDpT, expt.) and 
the scatter in the final PNC values is analyzed. The re- 
sults are summarized in Table I XII Our final value (cor- 
responding to data in Table |X]) is listed in the last row 
of Table IXII We note that essentially entire difference in 
the results comes from the dominant term (first row of 
Table m, and the variation in all other terms is insignif- 
icant. Therefore, the possible uncertainty in the next 
term (-1.35), which is bound to be substantial, can not 
be evaluated by this approach. While we have included 
the values with SDpT dipole matrix elements, there is 
no reason to expect these data to be more accurate than 
SD values. This conclusion is based on the breakdown of 
the correlation correction contributions and comparisons 
of the similar calculations in other monovalent systems 
that demonstrate cancellation of some missing effects in 
SD approximation but not in SDpT one. As a result, 
we conclude that the uncertainty in the dominant term 
owing to the Coulomb correlation correction is probably 
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TABLE XI: 'Scatter" analysis of the main part of the PNC 
amplitude (n = 7 — 10) in ^^^Ra"*". Lowest-order DP and 
random-phase RPA values are listed for reference. SD la- 
bels single-double all-order values, SDpT values include par- 
tial triple contributions. 



Energies 


{Mm 


{i\HpNc\j) 


TT^main 


DP 


DP 


DP 


38.95 


DP 


RPA 


RPA 


37.10 


SD 


SD 


SD 


39.05 


Expt. 


SD 


SD 


39.65 


Expt. 


SDpT 


SD 


40.22 


Expt. 


SD 


SDpT 


38.09 


Expt. 


SDpT 


SDpT 


38.65 


Expt. 


SD 


sd:. 


38.66 



" Scaled values are used for the 7pi/2 ~ 7s and 8pi/2 — 7s 
matrix elements only, remaining data are taken to be SD. 



TABLE XII: Contribution to the Epnc in ^^^Ra+ and com- 
parison with other theory. Our value for ^^^Ra^ is obtained 
by reducing our 223 value by 0.2% owing to the correc- 
tion for the different nuclear parameters and multiplying by 
138/135 neutron number ratio. All results are in units of 
W-''i\e\ao{-Qw/N). 



Isotope 


Term 


Value 


223 
223 
223 
223 
223 


TT-main 

TTitail 
-C'PNC 
TTiauto 

Breit 
Total 


38.66 
-0.02 
6.83 
-0.29 
45.18 


226 


Total 


46.09 


226 
226 


Mixed states [20] 
Sum over states [20] 
CCSD [21] 
CCSD(T) [21J 


42.9 
45.9 
46.1 
46.4 



on the order of 2%. We note that completely ab initio 
SD value is in good agreement with our final value. Mea- 
surement of the 6c?3/2 — 7pi/2 oscillator strength would 
help to reduce this uncertainty. 

We calculate remaining terms Ep^^c and i;|,^*g in both 
DF and RPA approximations. The RPA results are 
listed in Table IXIII together with our total value for 
the PNC amplitude. The corresponding DF results are 
£;auto ^ 4 8 ^tmi^ ^ The correction due to 

Breit interaction is obtained in the DF approximation. 
Our final value is compared with other calculations from 
Refs. [20l . [2l|. Our result for the terms with n < 7 and 
n > 9 (6.8) is in reasonably good agreement with the 
value from Ref. ^ (7.5). The notable feature of Ta- 
ble IXIII is an excellent agreement of all rather different 
high-precision calculations despite relatively large pos- 
sible uncertainties in various terms, with the exception 
of the mixed-states result [1^. Further calculations as 
well as experimental measurements will be necessary to 
achieve 1% accuracy in the PNC amplitude. 



TABLE XIII: Dependence of the lowest-order Ra+ PNC 
amplitude on the parameters of the nuclear distributions 
Cchargo(fm) and cpNc(fm). The parameter Cchargc is used in 
the charge distribution in the all-order wave function calcula- 
tions. The parameter cpNC is used in the modeling the nuclear 
density function in the PNC Hamiltonian. The variation of 
the given parameter is listed in % for convenience. The units 
for the PNC amplitude is 10""j|e|ao(-Qw/A^). 



^charge 


^L-cnargc 


CPNC 




^PNC 


"-tipNC 


D.oDl / 


no/ 
U/o 


D.oDl i 




44.y io 




6.8960 


0.5% 


6.8617 




44.853 


-0.13% 


6.9303 


1% 


6.8617 




44.792 


-0.27% 


6.9989 


2% 


6.8617 




44.671 


-0.54% 


7.2048 


5% 


6.8617 




44.310 


-1.34% 


6.8617 




6.8960 


0.5% 


44.875 


-0.08% 


6.8617 




6.9303 


1% 


44.837 


-0.17% 


6.8617 




6.9989 


2% 


44.761 


-0.34% 


6.8617 




7.2048 


5% 


44.531 


-0.85% 


6.8960 


0.5% 


6.8960 


0.5% 


44.815 


-0.22% 


6.9303 


1% 


6.9303 


1% 


44.717 


-0.44% 


6.9989 


2% 


6.9989 


2% 


44.523 


-0.87% 


7.2048 


5% 


7.2048 


5% 


43.954 


-2.14% 



We also investigated the dependence of the PNC ampli- 
tude on the values of the nuclear distribution parameters 
Cchargc and cpNC- As we described in the beginning of this 
section, the parameter Cchargc is used in the charge distri- 
bution in the all-order wave function calculations. The 
parameter cpNc is used in the modeling of the nuclear 
density function p(r) in the PNC Hamiltonian given by 
Eq. Both are modeled by the Fermi distributions; 

the all-order calculation is carried with both half-density 
parameters being equal to 6.8617 fm [sH- Since the DF 
result is rather close to the final value owing to various 
cancellations, it is sufficient to carry out this study us- 
ing DF data. The results are summarized in Table IXIII| 
where we list i?pNc calculated with varying values or ei- 
ther one or both parameters. The variation of the given 
parameter is listed in % for convenience. The results 
show that possible uncertainty in the PNC amplitude 
owing to the uncertainty in the value of the charge ra- 
dius (that is unlikely to be large) is negligible in compar- 
ison with the uncertainty in the correlation correction. 
For example, difference in the rms radii for A=223 and 
A=226 isotopes corresponds to the change in Cchargc that 
is on the order of 0.5% resulting in only 0.2% change 
in the PNC amplitude. Possible variance in the density 
p(r) in Eq. (|13p which is approximately neutron density is 
higher, but even 5% change in cpNC with the fixed value 
of the Cchargc Icads to 0.85% change in the PNC ampli- 
tude value. Table IXIIII may be used to recalculate the 
values of the PNC amplitude between different isotopes 
since the change in i?pNC with the nuclear parameters is 
essentially linear. 
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V. CONCLUSION 

We have calculated the energies, transition matrix el- 
ements, lifetimes, hyperfine constants, quadrupole mo- 
ments of the 6d states, as well as dipole and quadrupole 
ground state polarizabilities, and PNC amplitude in 
^^•^Ra^ using high-precision all-order method. The en- 
ergies of the 9pi/2, lOpi/2, and IOP3/2 levels are pre- 
dicted. The results for atomic properties are compared 
with available theoretical and experimental data. The 

6c?3/2 transition is found to 



PNC amplitude for the 7s - 



be 45.2 X 10-"i|e|ao(-Qiv/A^)- The dependence of the 
PNC amplitude on the choice of nuclear parameters is 
studied. This work provides a number of recommended 
values for yet unmeasured properties of Ra'^. 
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